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The layered compound a-RuCla is composed of a honeycomb lattice of magnetic Ru®'*' ions with 
the 4d® electronic state. We have investigated the magnetic properties of a-RuCla via magnetization 
and specific heat measurements using single crystals. It was observed that a-RuCla undergoes a 
structural phase transition at Tt ~ 150 K accompanied by fairly large hysteresis. This structural 
phase transition is expected to be similar to that observed in closely related CrCls. The magnetiza¬ 
tions and magnetic susceptibilities are strongly anisotropic, which mainly arise from the anisotropic 
g-factors, i.e., gab —2.5 and gfc —0.4 for magnetic fields parallel and perpendicular to the ab plane, 
respectively. These g-factors and the obtained entropy indicate that the effective spin of Ru®'*' is 
one-half, which results from the low-spin state. Specific heat data show that magnetic ordering 
occurs in four steps at zero magnetic field. The successive magnetic phase transitions should be 
ascribed to the competition among exchange interactions. The magnetic phase diagram for H || ab 
is obtained. We discuss the strongly anisotropic p-factors in a-RuCla and deduce that the exchange 
interaction is strongly XY-like. a-RuCla is magnetically described as a three-dimensionally coupled 
XY-like frustrated magnet on a honeycomb lattice. 

PACS numbers: 75.10.Jm, 75.30.Kz, 75.40.Cx 


I. INTRODUCTION 

It is known that a honeycomb-lattice antiferromag- 
net with the nearest-neighbor exchange interaction un¬ 
dergoes a conventional magnetic ordering even for the 
spin-1/2 case. However, when a certain amount of 
second- and third-neighbor exchange interactions or a 
certain amount of anisotropic exchange interaction ex¬ 
ists, the honeycomb-lattice quantum magnet exhibits an 
unusual ground state. In the last decade, spin-1/2 quan¬ 
tum magnets on honeycomb lattices have been attract¬ 
ing considerable attention from the viewpoints of the 
frustrated Ji — J 2 model— and the Kitaev-Heisenberg 
modelfi^— both of which can exhibit the spin liquid state 
in some parameter range. a-RuCb appears to be a spin- 
1/2 honeycomb-lattice magnet— Recently, great in¬ 
terest has been shown in the magnetic properties of a- 
RnCh^-^ 

a-RuCla has a layered structure. The crystal structure 
was first reported to be trigonal, P3il2^d^ but later it 
was found to be monoclinic, C2(m^ which is the same 
as the room-temperature crystal structure of CrClsj^i^ 
Figures [Ija) and (b) show the crystal structure of a- 
RuCls. The crystal structure is composed of RuCle oc- 
tahedra, which are linked in the ab plane by sharing 
edges. Magnetic Ru^+ ions with the icf’ electronic state 
form a slightly distorted honeycomb lattice. It has been 
reported that a-RuCla undergoes magnetic ordering at 
Tn = 13 — 15.6 K— However, little is known about 
the magnetic properties of a-RuCla. 

When Ru^’*' has the high-spin state, the total angular 
momentum and total spin are given by L = 0 and 5" = 5/2, 


respectively; thus, the magnetic moment is given by the 
spin only. Consequently, the exchange interaction and g- 
factor become isotropic. On the other hand, when Ru^+ 
has the low-spin state, the orbital and spin states are 
described by ^ = 1 and S'= 1/2, respectively. In this case, 
the magnetic moment is given by the effective spin- 1 / 2 , 
which is composed of the orbital angular momentum and 
true spin. In general, the exchange interaction and g- 
factor for the effective spin are fairly anisotropic. 

In this paper, we present the results of magnetization 
and specific heat measurements on a-RuCla. It was ob¬ 
served that the magnetic susceptibilities in a-RuCla are 
strongly anisotropic, which indicates the low-spin state 
of Ru^’*' together with small entropy. It was found that 
the magnetic ordering occurs in multiple steps, which 
we consider to be due to competition between exchange 
interactions. This paper is organized as follows. Experi¬ 
mental procedure is described in Section H. Experimental 
results are given in Section HI. The exchange interaction, 
the 5 -factor for the low-spin state of Ru^"*" and the phase 
diagram are discussed in Section IV. Section V is devoted 
to a conclusion. 


II. EXPERIMENTAL DETAILS 

Single crystals of RuCla were grown from a melt by the 
vertical Bridgman technique. Fine-grained RuCla was 
dehydrated in a quartz tube at 100 °C for three days. 
The temperature of the center of the furnace was set at 
1100 °C and the quartz tube was moved downward in 
the furnace at a rate of 3mm/h over 80 h. The crystals 
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FIG. 1: Monoclinic crystal structure {C2lm) of a-RuCla at 
room temperature (a) viewed perpendicular to the ab plane 
and (b) viewed along the b axis. Small purple and large green 
spheres are Ru^"*" and Cl“ ions, respectively. Thin solid lines 
denote the chemical unit cells, (c) Trigonal crystal structure 
(i?3) of CrCla in the low-temperature phase viewed along the c 
axis, which is expected to be the same as the low-temperature 
crystal structure of a-RuCla. 
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FIG. 2: Temperature dependence of (a) magnetic suscepti¬ 
bilities x~^/H and (b) inverse magnetic susceptibilities in 
a-RuCla measured at H = 0.1 T for H || ab and H T ab. The 
arrow indicates the structural phase transition temperature 
Ft. 


in magnetic fields of up to 9 T using a physical property 
measurement system (Quantum Design PPMS) by the 
relaxation method. 


III. EXPERIMENTAL RESULTS 



obtained were black and had wide surfaces parallel to the 
ab plane. The crystals were soft and easily bent like foil. 

Magnetization measurements were performed using a 
SQUID magnetometer (Quantum Design MPMS XL) in 
the temperature range 1.8 K < T < 100 K in magnetic 
fields of up to 7 T. Magnetic fields were applied parallel 
and perpendicular to the ab plane. High-field magne¬ 
tization measurement in a magnetic field of up to 57.5 
T was performed at 4.2 and 1.3 K using an induction 
method with a multilayer pulse magnet at the Institute 
for Solid State Physics, University of Tokyo. The abso¬ 
lute value of the high-field magnetization was calibrated 
with the magnetization measured by the SQUID magne¬ 
tometer. The specific heat was measured down to 0.36 K 


Figures[2Ka) and (b) show the temperature dependence 
of the magnetic susceptibilities and inverse susceptibili¬ 
ties measured for a magnetic field H parallel and per¬ 
pendicular to the ab plane. The magnetic susceptibil¬ 
ity for H \\ ab plane is much larger than that for H J-ab 
plane. The strongly anisotropic susceptibility is mainly 
attributed to the strongly anisotropic g-factor, which re¬ 
sults from the low-spin state of Ru^+. The anisotropy of 
the g factor is discussed in detail in the next section. 

In the inverse susceptibility for H J-ab plane, a discon¬ 
tinuous change is observed at Tt ~ 150 K. Figure [3] shows 
an enlargement of the magnetic susceptibility for H J-ab 
plane around Tf Hysteresis is clearly observed at the 
transition temperatures, Tt = 141 and 167 K upon cool- 
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FIG. 3: Hysteresis in magnetic susceptibility observed around 
r=150 K for H±ab. 




T[K] 

FIG. 4: (a) Low-temperature magnetic susceptibilities in q- 
RuGla measured at H^O.IT for (a) H || ab and (b) H _Lafe. 
Arrows indicate the anomalies caused by magnetic phase tran¬ 
sitions. 


ing and heating, respectively. This anomaly in the sus¬ 
ceptibility is ascribed to the structural phase transition. 
Although we performed X-ray crystal analysis below Tt, 
we could not determine the low-temperature structure 
because sharp X-ray spots were not observed owing to 
the softness of the crystal. In the closely related com¬ 
pound CrCla, a structural phase transition from the mon- 
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FIG. 5: Low-temperature magnetic susceptibilities in a- 
RuGls measured at various magnetic fields for H || ab. The 
susceptibility data are shifted upward by multiples of 2 x 10“® 
emu/mol. Arrows indicate the transition temperatures deter¬ 
mined from the specific heat anomaly, (b) Derivative of x with 
respect to T measured at 7L = 6 T. Arrows indicate the transi¬ 
tion temperatures determined from the specific heat anomaly, 
which are close to those giving local maxima or minima of 
dxjdT. 


oclinic structure {C2/m) to the trigonal structure (A3) 
takes place at Tt ~ 240 K.— Thus, it is likely that the low- 
temperature structure is the same as the low-temperature 
structure of CrCla, which is shown in Fig.[TKc). It is noted 
that the specific heat data show no sharp anomaly at the 
structural phase transition temperature Tt ~ 150 K. We 
consider that this is because the relaxation method used 
in this work is less sensitive to the first order phase tran¬ 
sition with large latent heat. 

Figure 0] shows the low-temperature magnetic sus¬ 
ceptibilities measured for H\\ab and H Tab. Clear 
anomalies indicative of magnetic phase transitions are 
observed at T = 13.9 and 7.6 K in both susceptibility 
data. These phase transitions are confirmed by perform¬ 
ing specific heat measurements as shown below. Our low- 
temperature magnetic susceptibilities are consistent with 
those reported recently by Sears et al.— and Majumder 
et alr^ However, we observed that for H T a5, the mag¬ 
netic susceptibilities in some a-RuCla samples decrease 
at 7.6 K with decreasing temperature in contrast to the 
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FIG. 6: Magnetic field dependence of the magnetization M 
and its field derivative dM/dH measured at 1.3 K for H || ab 
in a-RuCla for the highest fields of (a) 57.5 T and (b) 14 T. 
Arrows indicate the transition fields. 


behavior shown in Fig.|4jb). 

Figure [5lja) shows the low-temperature magnetic sus¬ 
ceptibility X = M /H measured at various magnetic fields 
for H II ab. Arrows indicate the transition temperatures 
determined from the specific heat anomaly. These transi¬ 
tion temperatures are close to those giving local maxima 
or minima of dx/dT, as shown in Fig. [5][b). The spike 
anomaly at 12 K is due to an instrumental problem. With 
increasing magnetic field, the transition temperatures de¬ 
crease, and for > 4 T, the high-temperature transition 
splits into two transitions. 

Figures [ 5 ] and [ 7 ] show the magnetic field dependence 
of magnetization and its field derivative dM/dH mea¬ 
sured at 1.3 K for H\\ah and H dLob, respectively. In 
both hgures, the upper and lower panels show data taken 
with the highest magnetic fields of 57.5 and 14 T, re¬ 
spectively. The absolute values of the magnetization for 
H \\ab and H dLob are considerably different, as observed 
in the magnetic susceptibilities. This is ascribed to the 
strongly anisotropic (/-factors. Extrapolating the magne¬ 
tization curves to higher fields, we estimate the (/-factors 
to be (/ab = 2.5±0.2 and (/c = 0.40 ± 0.03. As shown by 
arrows in Figs. [ 5 ] and [71 some anomalies in dM/dH in¬ 
dicative of field-induced phase transitions are observed at 


FIG. 7: Magnetic field dependence of the magnetization M 
and its field derivative dM/dH measured at 1.3 K for H J^ab 
in (r-RuCls with the highest fields of (a) 57.5 T and (b) 20 T. 
Arrows indicate the transition fields. 


= 1.6, 6.2, 7.5, 7.9 and 9.6 T for H || ab, and at 77 = 7.9 
and 9.1 T for 77 T ab. For 77 || ab, the field-induced phase 
transition occurs in many steps. However, a distinct tran¬ 
sition to saturation is not observed for either field direc¬ 
tion, despite the sufficiently low temperature of 1.3 K. 
This indicates that the total spin is not conserved. Be¬ 
cause a-RuCla is considered to be a localized spin system 
at helium temperatures, a strong antisymmetric interac¬ 
tion such as the Dzyaloshinskii-Moriya (DM) interaction 
may be responsible for the smearing of the saturation 
transition. 

Figure [5] shows the temperature dependence of specific 
heat C and C/T measured at various magnetic fields for 
77 II ab. The anomaly in C/T near 14 K is due to the 
instrumental problem. Our specific heat data are con¬ 
sistent with those reported by Majumder et ah^. At 
zero magnetic field, a sharp peak is observed at rN 4 = 7.6 
K. Above 7.6 K, three small anomalies are observed at 
Tns = 10.4 K, Tn 2 = 12.3 K and Tni = 13.8 K as indicated 
by arrows. We consider that these anomalies arise from 
magnetic phase transitions, because they are sharper 
than those owing to the short-range spin correlation. The 
present result shows that the magnetic ordering occurs 
in four steps. The transitions at Tni and Tn 4 are also 
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FIG. 8: Low-temperature specific heat C and C/T in a-RuCls measured at various magnetic fields for H || ab. The specific 
heat data and C/T are shifted upward by multiples of 2 J/mol-K and 0.2 J/mol-K^, respectively. Arrows show the anomalies 
indicative of magnetic phase transitions. 


clearly observed in the magnetic susceptibility as shown 
in Fig. m 



FIG. 9: Temperature dependence of Cmag/T in a-RuGla mea¬ 
sured at zero magnetic field. The solid curve represents the 
magnetic entropy Smag- 

The magnetic specific heat Cmag was evaluated using 
the specific heat of nonmagnetic ScCla, which has a sim¬ 
ilar layered crystal structure to RuCla.— Figure [9] shows 
the temperature dependence of Cmag/T and magnetic en¬ 
tropy S'mag at zero magnetic field. C^aglT has a broad 
maximum around 85 K, which is interpreted to be caused 
by the short-range spin correlation. The magnetic en¬ 
tropy obtained at 140 K (<Tt) is S’mag = 3.8 J/mol-K, 
which is approximately two-third of i? In 2 = 5.76. This 
small magnetic entropy is consistent with the low-spin 
state of Ru^+ with effective spin-I/2. 

As shown in Fig. [51 with increasing magnetic field 
for H \\ ab, the peak heights of the anomalous specific 
heat for all the phase transitions decrease and all the 
transition temperatures shift gradually towards the low- 
temperature side. The lowest transition temperature Tn 4 



H(T) 


FIG. 10: Field scans of the specific heat in o-RuCla measured 
at 1.3, 2.4 and 3.5 K for H || ab. Arrows show the anomalies 
indicative of magnetic phase transitions. 


decreases rapidly above 7 T. This can be observed in the 
field scans of the specific heat at three temperatures be¬ 
low 4 K, as shown in Fig. (TUI The peak around 7.8 T, 
indicative of a magnetic phase transition, is interpreted 
to be connected to the transition line for Tn 4 for H <7 
T. The high-temperature transition T^i splits into two 
transitions above H = 5 T. 

The transition data obtained for H \\ ab are summa¬ 
rized in Fig. |TT] Four ordered phases exist at zero mag¬ 
netic field and six ground states exist in magnetic fields. 
The features of a-RuClg observed in the present measure¬ 
ments are strongly anisotropic magnetic properties and 
many ordered phases. In the next section, we discuss 
these features. 
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FIG. 11: Magnetic phase diagram for H || ab. The circular 
and triangular symbols are transition points determined from 
the temperature and field dependences of specific heat, re¬ 
spectively, and the rectangular symbols are those determined 
from the magnetization process. Solid and dashed lines are 
visual guides. 


IV. DISCUSSION 


A. Effective exchange model and p-factor 


In this subsection, we discuss the strongly anisotropic 
magnetic properties and derive an effective model that 
describes the low-temperature and low-energy magnetic 
properties of a-RuCla in accordance with the effective 
model of Co^+ in an octahedral environmentBe¬ 
cause the magnetic entropy obtained from the specific 
heat data indicates that the effective spin of is one- 
half, it is natural to assume that the five electrons in the 
Ad orbitals of Ru^+ are in the low-spin state owing to the 
strong crystalline helds from surrounding Cl” ions. 

In the low-spin state, all five electrons in the Ad orbitals 
occupy the de orbital. Because the matrix elements of 
the orbital angular momenta and with respect 

to the orbital states and (jjQ for the de orbital are 

given by changing the sign of those for the orbital an¬ 
gular momenta and with respect to the orbital 

states (j)x,(t>v s-nd 4>z for the p orbital, respectively, we 
can replace Ide by —I with 1 = 1. The spin-orbit coupling 
of these electrons is expressed as 


^SO — 2_^ 3 ^ 3 {li-Si). ( 1 ) 

i—1 ^ ^ 


For the three electrons with up spin, their orbital angu¬ 
lar momenta cancel out, = ih = 0- For the other two 
electrons with down spin, their spin {i = A and 5) is ex¬ 
pressed using the total spin S with 5= 1/2 as = — S. 
For these reasons, the spin-orbit coupling of eq. is 


written as 

nso = gdB^z(^^'^{l-S) = X{l-S). ( 2 ) 

The coupling constant A is positive and its magnitude has 
been reported to be A ~ 1000 cm”^i^ When the p orbital 
of the surrounding Cl” is mixed with the Ad orbitals of 
Ru^“*", the matrix elements of the angular momentum I 
are reduced. This effect is expressed by replacing I with 
kl with 0 < fc < 1. 

The orbital state of the low-spin state of Ru^"*" in an 
octahedral environment is triply degenerate. The orbital 
degeneracy can be lifted by the spin-orbit coupling and 
the trigonal crystalline field, which are written as 

n' = X'{l-S) + S{{Ff-2/3}, (3) 

where A' = kX, and the second term represents the en¬ 
ergy of the trigonal crystalline field. When the RuCle 
octahedron is trigonally compressed, i5 > 0, and when it 
is elongated, <5 < 0. The orbital triplet splits into three 
Kramers doublets. Their eigenvalues are expressed as 

A' 3A' 2’ ^ ’ 

and 



These eigenvalues are shown in Fig.[T2] as a function of 
(5/A'. 

When the temperature T is much lower than 
A'~ 1000 cm”^^ i.e., T< lOOK, the magnetic prop¬ 
erty is determined by the lowest Kramers doublet with 
E = E~. The eigen-states of the lowest Kramers doublet 
are expressed as 

V-± =ci|±1,t1/2) + C2|0,±1/2), (6) 



FIG. 12: Energy levels of Ei and as a function of 5/X'. 
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where \mi,ms) denotes the state with l^=mi and When a magnetic field is applied parallel to the trigonal 
=ms- Coefficients Ci and C 2 are given by axis, the Zeeman energy is expressed as 


Cl = ^1/1- 


V2y 


C 2 —- 7^ 1 1 + 


A 


V2V yA 2 +T 


with 


A = 


2(S/A') - 1 
2^2 


, ( 7 ) 


( 8 ) 


Within the lowest Kramers doublet, we have 
(V'il'S'liAi) = T^cf - ci), 
{^|J+\S+\^|;-) = {y;-\S-\yj+) = 4 . 


Using these relations, we can replace the true spin S with 
5"= 1/2 by the spin- 1/2 operator s given by 


nl = (14) 

with 

= 2\{{k + 1)4 - cl}\. (15) 

Note that the Zeeman levels of ip- and ip+ are reversed 

at (c 2 /ci)^ = fc-|-1. The Zeeman energy for a magnetic 
field perpendicular to the trigonal axis is expressed as 

+1~) + 2(^+ + S-)\iP^)H 
= (16) 


with 


s^ = 4 s 4 sy = clsy, =-(4-4)3^. (lo) 

We assume that the exchange interaction between true 
spins Si and Sj is described by the Heisenberg model 
T^ex = JSi-Sj. Substituting eq. (flUl) into Hex, we obtain 
the effective model 

HeS = (S/Sj + s/sj) + (11) 

with 

= 4j, jii = (4 - 4)^j. ( 12 ) 

The exchange constants J-^ and in the special cases 
are shown in Table HI When the trigonal crystalline field 
is absent (d = 0 ), the effective exchange interaction HeS 
becomes the Heisenberg model, while when (5/A'= 1/2, 
it becomes the complete XY model. For i5/A'—>- 00 , 
the orbital angular momentum is quenched, so that the 
magnetic moment is given by the spin only, and HeS 
again becomes the Heisenberg model. For 5/X' < 0, HeS 
is Ising-like and becomes the complete Ising model for 
(5/A' —>■ — 00 . 


TABLE I: Coefficients ci and C2, exchange constants and 
and p-factors and in the special cases of (5/A'. 


S/X' 

Cl 

C2 


Jll 

9^ 

7' 

0 


-1/73 

J/9 

J/9 

2(l + 2fc)/3 

2(1 -h2fc)/3 

1/2 

1/72 

-1/72 

J/4 

0 

1 -t 72fc 

k 

1 

1/73 

4J/9 

J/9 

4(l-Hfc)/3 

2(1 - fc)/3 

00 

0 

-1 

J 

J 

2 

2 

— 00 

1 

0 

0 

J 

0 

2{k + l) 


The lowest Kramers doublet splits into two Zeeman 
levels when subjected to a magnetic field. The splitting 
of the Zeeman levels in Ru^+ has been discussed by many 
authorsi^^— The Zeeman term is written as 

(13) 


g^ = 2(4 - V2kciC2). ( 17 ) 

Figure [ 13 ] shows these ^-factors as a function of (5/A'. The 
(/-factors in the special cases are shown in Table HI 

The (/-factors estimated from the high-field magnetiza¬ 
tion process are gab = g^p — 2.5 and gc = glxp — 0.4. This 
indicates that (5/A' ^ 1 in RuCla. Figure [TST bl shows the 
behavior of the (/-factors in the range of 0 < 6/X' < 2 , 
where g^^ changes rapidly with varying S/X'. There 
are two sets of parameters, ((5/A', fc) = (0.77, 0.95) and 
(1.18, 0.83), that satisfy the (/-factors observed for a- 
RuCls. In the present experiment, it is difficult to evalu¬ 
ate which set of parameters is realized for RuCls. The ex¬ 
change anisotropy is given as = 0.099 and 0.37 for 

(5/A'= 0.77 and 1.18, respectively. Hence, the exchange 
interaction between effective spins is strongly XY-like in 
cr-RuCla. From the above discussion, we can conclude 
that the strongly anisotropic magnetic properties ob¬ 
served in a-RuCls arise from the trigonal crystalline field, 
which is close to the spin-orbit coupling, S/X' ~ 1. When 
uniaxial pressure is applied parallel to the trigonal axis, 
that is, normal to the ab plane, the coefficient S of the 
trigonal crystalline field will increase because the mag¬ 
nitude of the trigonal compression of the RuCle octahe¬ 
dron is increased by the uniaxial pressure. If S/X' = 0.77 
at ambient pressure, gc will decrease under the uniax¬ 
ial pressure, while if (5/A'= 1.18, gc will increase. Thus, 
magnetization measurements under the uniaxial pressure 
will be useful in determining the parameter S/X'. 

B. Successive magnetic phase transitions 

As shown in the phase diagram of Fig.[TTl cr-RuCla un¬ 
dergoes four magnetic phase transitions at zero magnetic 
field. If the frustration is absent, then a single magnetic 
phase transition will occur because the honeycomb lattice 
is bipartite. Hence, it is considered that the successive 
phase transitions arise from the frustration effect owing 
to the competing interactions. 


Hz = -^lB{-kl + 2S)■H. 
















FIG. 13: (a) g-factors as a function of 5/\'. (b) Enlargement 
of the g-factors between 5/A' = 0 and 2. The two horizontal 
lines are experimental g-factors estimated from the high-held 
magnetization process. Open circles and squares represent 
two sets of the factors suitable for a-RuCls. 

Recently, Sears et alr^ investigated the spin struc¬ 
ture below Tisi 4 = 7.6 K at zero magnetic field by neutron 
diffraction. They reported that the so-called zigzag-type 
order is realized below Tn 4 . They discussed the ground 
state from the viewpoint of the competition between the 
Heisenberg term JSi ■ Sj and the Kitaev term —KS^Sj, 
where 7 corresponds to the direction of the bond con¬ 
necting Si and Sj. A possible origin of the Kitaev term 
is discussed in Appendix. However, it appears difficult to 
derive the successive phase transitions within the nearest 
neighbor interactions. 

In MnBr 2 , Mnl 2 and NiBr 2 , which have the sim¬ 
ilar layered crystal structure to a-RuCla, two mag¬ 
netic phase transitions have been observed at zero mag¬ 
netic fieldi^^— Their successive phase transitions have 
been explained theoretically in terms of the competition 
among the exchange interactions up to the third neighbor 
and the interlayer exchange interactioni^^iM 

When MXe octahedra centered by magnetic ions 
M form a honeycomb or triangular lattice by sharing 
their edges, the nearest-neighbor exchange bond angle 
M—X—M is close to 90°. In a-RuCla, the nearest- 


neighbor exchange bond angle Ru—Cl—Ru is approxi¬ 
mately 96°. When the bond angle is close to 90°, the 
exchange tends to be ferromagnetic or weak, even if it is 
antiferromagnetic 1 ^ Actually, in closely related CrCls, 
a ferromagnetic ordering is realized in the honeycomb 
latticej^i^ The second- and third-neighbor exchange in¬ 
teractions J 2 and J 3 in a-RuCla are considered to be 
the same order of magnitude as the nearest-neighbor ex¬ 
change interaction Ji, as observed in MnBr 2 , Mnl 2 and 
NiBr2. 

The classical ground state for the Ji — J 2 — J 3 Heisen¬ 
berg and XY model was theoretically investigated 
by Rastelli et al^ and Fouet et al^ They showed 
that zigzag ordering emerges when J 2 lJi>ll 2 and 
>7^3 /Ji>1/2 for antiferromagnetic Ji (> 0) and when 
« 72 /</i< 1/2 and J 3 /Ji <0 for ferromagnetic Ji(< 0 ). 
Thus, zigzag ordering is possible in a realistic parameter 
range. Therefore, we infer that the successive phase tran¬ 
sitions and zigzag magnetic ordering observed in a-RuCla 
are attributed to the competition among the exchange 
interactions up to the third neighbor and the interlayer 
exchange interaction. However, a theoretical description 
of the successive phase transitions is an open problem. 


V. CONCLUSION 

We have presented the results of magnetization and 
specific heat measurements on the honeycomb-lattice 
magnet a-RuCla. This compound undergoes a first-order 
structural phase transition at Tt = 154 ± 13 K. The struc¬ 
tural phase transition is expected to be a transition from 
the monoclinic room-temperature structure ((72/771) to 
a trigonal structure by analogy with that observed in 
closely related CrClai^ 

The magnetic susceptibility and magnetization are 
strongly anisotropic, i.e., these quantities for H || ab are 
much larger than those for H T ab. This is ascribed 
to the strongly anisotropic g-factor characteristic of the 
low-spin state of Ru^"*" with the electronic state. 
We discussed the effective exchange interaction and g- 
factor taking the spin-orbit coupling and trigonal crys¬ 
talline field into consideration. We demonstrated that 
the strongly anisotropic magnetic properties observed in 
a-RuCla occur when the magnitudes of the spin-orbit 
coupling and trigonal crystalline field are close to each 
other, i.e., 5/A'~l, and that the effective exchange in¬ 
teraction is strongly XY-like in a-RuCla. 

It was found from the specific heat and magnetization 
measurements that a-RuCla undergoes four magnetic 
phase transitions at zero magnetic field and five field- 
induced transitions at T = 0. We presented the magnetic 
field vs temperature phase diagram for H || ab in Fig. 1111 
We suggest that these successive phase transitions are at¬ 
tributed to the competition among the nearest-neighbor, 
second- and third-neighbor exchange interactions. 
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FIG. A14: Classical configuration of spins Si and Sj and co¬ 
ordinate systems. The y direction is chosen to be parallel to 
the bond vector ri 2 and D. 
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Appendix A: Possible origin of the Kitaev term 

Here, we give a possible origin of the Kitaev term. As 
discussed above, the dominant exchange interaction be¬ 
tween magnetic moments is the spin-1/2 XXZ model ex¬ 
pressed by eq. dm, which is isotropic in the ah plane. 
We consider the Dzyaloshinskii-Moriya (DM) interaction 
D ■ [si XS 2 ] with the D vector parallel to the bond vector 
't '12 =i '2 — ri. This condition is allowed when the space 
group is P3il2. In this case, there is a twofold axis pass¬ 
ing two neighboring spins, which leads to D || ri 2 .— The 
interaction between the nearest-neighbor spins si and S 2 
is expressed as 

ni2 = J^{slsl + s\sl) + jhls^^+D-[sixS2], (Al) 


where the z axis is chosen to be normal to the honey¬ 
comb lattice. When the XXZ-type exchange interaction 
is antiferromagnetic, the stable classical spin configura¬ 
tion is a canted antiferromagnetic state as illustrated in 
Fig. IA14I The canting angle 0 is given by 

2D 

= (A2) 

Here, we define the local coordinates x'yz'—O and 
x"yz"—0 as shown in Fig. IA141 i.e., the x' and x" axes 
are taken to be parallel and antiparallel to the spins Si 
and S 2 , respectively. The y axis is parallel to the bond 
vector and D. The spin operators 5^2 sf 2 in the 
original coordinate system are expressed as 

Si = sf cos 6 — sf sind, sf = sf sin 6 + sf cos 0 , 
sf = sf cos 0 ^ sf sin 0, S 2 = —S 2 sin 0 + sf cos 0. 

(A3) 

Substituting eq. (IA3II into eq. dni, we obtain 

Hi 2 = {J^ + K)(sfsf + s\ si) + (Jll + K)sfsf 

- Ks\sl (A4) 

with 

K = 2D sin 0 cos 0 — {J^ + ) sin^ 0. (A5) 

Because the y axis is parallel to the bond vector ri 2 , 
which has three directions in the crystal depending on 
the configuration of Si and S 2 , eq. (IA4I) is equivalent to 
the spin-1/2 Kitaev-XXZ model. As described above, the 
DM interaction can be the origin of the Kitaev term. 
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